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Abstract 

We study Hermitian unitary matrices S £ C"'" with the following property: There exist 
r > and t > such that the entries of S satisfy \Sjj \ = r and \Sjk \ ~ t for all j, fc = 1, . . . , n, 
j 7^ k. We derive necessary conditions on the ratio d :— r/t and show that they are very 
restrictive except for the case when n is even and the sum of the diagonal elements of S is 
zero. Examples of families of matrices S are constructed for d belonging to certain intervals. 
The case of real matrices <S is examined in more detail. It is demonstrated that a real S can 
exist only for d = — 1, or for even n's and d satisfying ^ + d = 1 (mod 2). We provide 
a detailed description of the structure of real S with d > f — |, and derive a sufficient and 
necessary condition of their existence in terms of the existence of certain symmetric {v, k, A)- 
designs. We prove that there exist no real S with d £ (f' ~ !> j ~ §)■ ^ parametrization of 
Hermitian unitary matrices is also proposed, and its generalization to general unitary matrices 
is given. At the end of the paper, the role of the studied matrices in quantum mechanics on 
graphs is briefly explained. 



1 Introduction 

Unitary matrices with various special properties emerge in a wide scale of applications in physics 
and in the engineering, and at the same time they constantly attract the attention of pure math- 
ematicians. One of the most fascinating and longest-standing problems in mathematics is the 
Hadamard conjecture: // n is a multiple of 4, then there exists an n x n matrix H with entries 
from {—1, 1} such that HH^ = nl . Although the conjecture is believed to be true, no proof has 
yet been found. The matrix H with these properties is called Hadamard matrix of order n, and 
is just a multiple of an orthogonal matrix having all the entries of the same moduli. Hadamard 
matrices have numerous practical applications in coding, cryptography, signal processing, artificial 
neural networks and many other fields, see e.g. the monography [1]. 

A similar problem is related to the existence of so-called conference matrices. A conference 
matrix of order n is an n x n matrix C with on the diagonal and ±1 off the diagonal such that 
CC^ = (n — 1)/. Matrices of this type are important for example in telephony and in statistics, 
but as in the case of Hadamard matrices, there is still no definite characterization of orders n for 
which a conference matrix exists. 

Note that both Hadamard and conference matrices have these two properties: 



(PI) they are multiples of orthogonal matrices; 

(P2) all their off-diagonal entries are of the same moduli, and also all their diagonal entries are of 
the same moduli. 

These properties can serve as an inspiration to generalize Hadamard and conference matrices to the 
whole set of matrices satisfying (PI) and (P2). A subclass fulfilling a certain additional condition, 
namely the class of matrices with constant diagonal, has been studied in [2, 3]. 

Both Hadamard and conference matrices are by definition real, but they can be naturally 
generalized to complex ones by allowing their entries to take any values from the unit circle instead 
of {1, —1}. Complex Hadamard and conference matrices and their properties are nowadays widely 
studied as well, see e.g. [1, 4]. This fact may serve as another inspiration for generalizations: 
Examine all unitary matrices satisfying (P2). 

The subject to be discussed in this paper is close to the aforementioned generalization. We 
will study complex unitary matrices satisfying (P2) that are at the same time Hermitian. Our 
aim is to examine their existence and their properties, and perhaps to motivate a more extensive 
study of them, as they play an important role in the quantum mechanics on graphs (we will devote 
Section 7 at the end of the paper to a more detailed explanation). Since the real matrices of this 
type are for many reasons interesting, we will focus on the real case in a separate section. Another 
purpose of the paper is to propose a parametrization of unitary matrices, with a particular accent 
put on their Hermitian subset. 



2 Preliminaries 

Definition 2.1. (i) A square matrix M € C"'" is called permutation- symmetric if there are 
a, 6 e C such that the entries of M satisfy 

Mjj = a and Mjk = b for all j, fc = 1, . . . ,n, j ^ k . 

(a) We call a square matrix M G C"'" modularly permutation-symmetric if there are a,b > 
such that the entries of M satisfy 

\Mjj\ = a and \Mjk\ = b for all j,k = 1, . . . ,n, j ^ k . 

"Modularly permutation-symmetric" will be hereinafter abbreviated as MPS. Matrices from 
Definition 2.1 have the following property: If M is a permutation-symmetric matrix (or an MPS 
matrix) and P is a permutation matrix of tlie same size, then PMP^^ is a permutation-symmetric 
matrix (or an MPS matrix, respectively) as well. 

In this paper we are particularly interested in unitary and at the same time Hermitian modu- 
larly permutation-symmetric matrices; we will denote them by the symbol S. As diagonal Hermi- 
tian unitary MPS matrices are trivially of the form S = diag(±l, ±1, . . . , ±1), from now on we will 
focus on the case when the modulus of the off-diagonal entries is nonzero. For the sake of brevity, 
let us denote the set of all Hermitian unitary MPS matrices with the ratio d := ^oe^df^gonaTentiyl 
by the symbol A4n{d), i.e., 



Mn{d) = \S& U(n) 



S is MPS A = d A S = S* 
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in other words, elements of A4n{d) are Hermitian unitary matrices n x n of the type 

giai2 e'"i^ • • • e'^i" \ 
±d e'"23 ... c'"=^" 

e-ia23 ±(l ... Qias^ 



Q-i0C2n Q-i0C3n ... ±d J 

Remark 2.2. For d = and d= 1, Mn{d) represents the set oinxn Hermitian conference matrices 
and Hermitian Hadamard matrices, respectively: 

• <S e A^Ti(O) iff C := \/rin • <S is a (complex) Hermitian conference matrix; 

• 5 e A^„(l) iff H := ^Jn • 5 is a (complex) Hermitian Hadamard matrix. 

Within each set A1„((i) wc introduce an equivalence: 

Definition 2.3. We say that matrices <Si,iS2 G M.n{d) are equivalent, written as <Si ~ <S2, if one 
can be obtained from the other by performing a finite sequence of the following operations: 

• for a certain j, fc, transpose the j-th and the fc-th row, and at the same time transpose the 
j'-th and the fc-th column; 

• for a certain j and G K, multiply the j-th row by e'"^, and at the same time multiply the 
j-th column by e~"^; 

• multiply the whole matrix by —1. 

In other words, S\ ^ Si iff there cixist a permutation matrix P and a diagonal unitary matrix 
D = diag(e''^i , e'"^- ) such that 

Si = DPSiP-^D-^ or Si = -DPS2P~^D-^ . 

Remark 2.4. In the literature on Hadamard matrices, a weaker equivalence is mostly used, namely 
that the operations can be performed independently on the rows and on the columns. Within the 
set M.n{d), however, we require the equivalence as it is defined above, mainly because it ensures 
the property (<Si e Mn{d) A<Si S2) <S2 e Mn{d). 

Notation 2.5. Everywhere in the paper, the symbols I^*^ and J*^*^) denote the identity matrix of 
order k and the matrix fc x fc all of whose entries are 1, respectively. 

Finally, let us give the definition of the symmetric {v, k, A)-design which will be useful for 
contructions of matrices <S G M.n{d) in Section 5 and at the end of Section 6. 

Definition 2.6. Let ?;>fc>A>lbe integers. A symmetric {v, k, X)-design is a pair V = {P, B), 

where P = {pi, . . . is a set of v points and B = {Bi, By} is a set of v subsets oiP (blocks) 
each containing k points, such that each pair of distinct points is contained in exactly A blocks. 

An incidence m,atrix A = {Aij) of P is a w x w matrix with entries from {0, 1}, where Aij = 1 
if and only if pj G B^. 

An A e {0, 1}^'" is an incidence matrix of a symmetric {v, k, A)-design if and only if 

AA^ = {k- X)!^"^ + AJ^") and AJ^") = AJ^") , (1) 

cf. [5], Thm. 2.8, or [1]. 



( ±d 
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3 Parametrization of unitary matrices 



This section addresses the problem of parametrization of unitary matrices. The result will be 
useful later in this paper, but we believe that it may be generally of interest in itself. The solution 
we present is based on ideas from [6] and [7]. 

We begin with the case when U G U(n) is Hermitian, and then we will generalize the parametriza- 
tion to all unitary matrices. At the end of the section it will be shown that after a certain minor 
upgrade, the parametrization is applicable much more generally, namely to Hermitian matrices H 
solving the equation H'^ = al + bH. 

Observation 3.1. If a matrix S is unitary and Hermitian, then the eigenvalues ofS are from the 
set {-1,1}. 

The most important result of this section follows. 

Theorem 3.2. (i) Let S be a Hermitian unitary matrix of order n. If S ^ ±7^"), then there 
exist an m £ {1, . . . ,n — 1} , a matrix T e C™'""™ and a permutation matrix P such that 

j{m) 



S = -7(") +2p( j (7('") + TT*) ^ ( 7('") T ) P"^ 



(2) 

_ ^ / -7(") +2 (7(") +Tr*)"^ 2 (7('") +Tr*)"^T \ p-^ 

\ 2T* (7('") + J^J^*'-j^^ _jin-m) _|_ ^/(™) 4^ J 

(a) For any to <E {1, . . . , n — 1}, for any T G £m,n-m permutation matrix P of order 

n, the matrix S given by (2) is Hermitian unitary. 

(Hi) If S is given by (2), then the columns of the matrices 

/("') \ { T 

rj.* 1 and P[ _j-(n-m) 

are eigenvectors of S corresponding to the eigenvalues 1 and —1, respectively. 

Proof, (i) Let 5 be a Hermitian unitary n x n matrix different from ±7'"^ and m denote the 
multiplicity of its eigenvalue 1. Since <S ^ ±7^"^ it holds to ^ and m^n, thus to S {1, ... , n—1}. 
The multiplicity of the eigenvalue —1 equals n — m, and therefore 

rank(5 + 7^")) =n- dimKer(5 + 7^")) = n - (n - m) = m e {1, • • • , n - 1} . 

Hence there is a regular M G C"'"* and a permutation matrix P such that 

^ \ T2M T2MT1 

note that P can be omitted if and only if the upper left submatrix to x m of <S + 7^"^ is regular. 
As 5 = 5*, necessarily A7 = M* and (MTi)* = T2M. Since M is regular, we have T2 = Tf. Let 
us set for brevity T :=Ti. Since <S is unitary, the matrix 

cc*-T{n),^pf M MT \p-i.p{ M^ + MTT*M M^T + MTT*MT > _ 

\T*M T*MT ) ^ \T*M'^ + T*MTT*M T* M^T + T* MTT* MT ' 

equals 7^"), hence we obtain 2M + + MTT*M = 0, equivalently 2M-^ = 7('") + TT*. Conse- 
quently, M = 2(7(™) + rr*)-i. 
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(a) Any iS given by (2) obviously satisfies SS* = /("^ and <S = <S*. 
(iii) If S is given by (2), a straightforward calculation gives 

SP ) = P ( ^^"^ ) and SP ( ) = ( -/("-™) 

therefore (iii) holds true. □ 

Remark 3.3. Let 5 ^ ±7^"^ be a Hermitian unitary matrix of order n, m = rank(»S + /^"^), and 
^(iii) be the upper left m x m submatrix of 5. It follows from the proof of Theorem 3.2 that the 
permutation matrix P must be involved in the parametrization (2) if and only if »S^^'^^ + I^"^^ is 
singular. 

Remark 3.4. Since the matrix T occurring in the parametrization (2) determines the eigenvectors 
of <S, it is related to the diagonalization of S as well. It follows from Theorem 3.2 (iii) that 



for 

Jim) J, N z - ( ^^"''^ ^ 
rp^ j{n—ra) I ' ^ I q jin—m) 



Xm. — P 



The main idea of Theorem 3.2 can be extended to a general unitary matrix: 

Theorem 3.5. Let U e U(n) such that U ^ — /(") . Let n — m denote the multiplicity of its 
eigenvalue —1. 

• If n — m ^ 0, then there exists a T G £,rn,n-m^ ^ Hermitian S G C"*'™ and a permutation 
matrix P such that 

jim) 




+TT* +\S)~^ (I^""^ + TT* + iSy^ T 
T* (7(™) + TT* +iS)~^ T* (7(™) + TT* +iS)~^T 



(3) 



and conversely, any matrix given by (3) is unitary. 

• //n - m = 0, there exists a Hermitian S e C'"'™ such that U = + 2 (J^™) + iS')~\ 

and conversely, any U given by this formula is unitary. 

Proof. Let n — 0. Similarly as in the proof of Theorem 3.2, we start from the decomposition 

[/+/(") = P M { I Ti )p-\ where M e C™'™ is regular, and then require UU* = /("). 

It leads to T2 — Tj* and M = 2(7 + TiT^ + iS)^^ for a certain Hermitian matrix S. 

If n — m = 0, the matrix U + J*^"-' is regular. Let us denote U + =: M. Then the requirement 

UU* = gives M = 2(7(") + 15)-^ for a certain Hermitian S. □ 

Remark 3.6. The idea from Remark 3.3 applies to (3) as well. The permutation matrix P must be 
involved in (3) iff U^^'^^ _(- is singular, where U^^'^^ stands for the upper left submatrix mxm 

of U and m = rank(U + /(")). In case U^^'^"^ + /*■'"•' is regular, P may be omitted. 
Remark 3.7. The imitary group U(n) has real parameters. There exist several known parame- 
trizations, i.e., ways how the parameters can be assigned to matrices U G U(n), for example [8, 9] 
and many other. In accordance with P. Dita (cf. e.g. [10]), we call a parametrization naiural if the 
involved parameters are free, i.e., there are no supplementary restrictions upon them to enforce 
unitarity. Our solution (3) falls within that class. On the other hand, (3) has a disadvantage that 
if the rows and columns of U are not suitably ordered, then a permutation matrix must be brought 
in, see Remark 3.6. 
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Hermitian solutions of quadratic matrix equations 



The reader may have observed in the proof of Theorem 3.2 that the essential properties of S 
that allowed us to obtain the parametrization (2) were the following two: the hermiticity of <S 
and the fact that <S has only two eigenvalues. In the light of this idea, we will generalize the 
parametrization (2), originally developped for Hermitian unitary matrices (i.e., solutions of <S^ = I), 
to Hermitian solutions of more general matrix quadratic equations 

H'^ = al + bH (a, 6 e K) . (4) 

We observe at first that the eigenvalues of any solution H of (4) must satisfy = a + bX, hence 
a{H) = {Ai,A2} where Ai,2 = ^ {b±Vb^ + 4a). Since A1.2 are real due to the hermiticity of H, 
one has to assume a,b E M. and 4a + b^ > 0. Note that the case 4a + &^ = is not interesting, 
because it represents the situation when any Hermitian solution of (4) has the eigenvalue b/2 with 
multiplicity n, thus H ^ ^I. For these reasons we shall assume the strict inequality 4a + 6^ > 0. 

Let us transform Equation (4) into its equivalent form 

by f b^ 



and define 

V4a + 6^ V 2 

Matrix M is Hermitian (because H is Hermitian) and at the same time unitary, since it satisfies 
= I . Therefore we can apply Theorem 3.2 and in this way obtain the sought parametriza- 
tion of i7, sec Theorem 3.8 below. We remark that the trivial solutions of (4), namely H = 
i (6 ± \/4a + b'^) I, are excluded from the parametrization, just as <S = ±/ have been excluded in 
Theorem 3.2. 

Theorem 3.8. Let a, & G R, 4a + 5^ > 0. 

(i) A Hermitian nx n matrix H different from \{b± v'4a + 6^) I*^") satisfies H"^ = al^^^ + bH 
if and only if 

H= ^~ ^^f+^ lin) + • P ( ) (J(-) + TT*)" ( J(-) T)P-' (5) 

for an m E {1, . . . , n — 1} , a m,atrix T G c™^"-™ ^ji^i q permutation matrix P. 

(a) If H is given by (5), then the columns of the matrices 

j(m) \ ( T 

J,* j and P\ _j(n-m) 

are eigenvectors of H corresponding to the eigenvalues ''+V4a+b^ ^^.^^ b-V4a+b^ ^ respectively. 



4 Modular permutation symmetry 

In the following part of the paper we will study Hermitian unitary MPS matrices. Prior to that, 
let us bring in a proposition characterizing the set of Hermitian unitary permutation- symmetric 
matrices (cf. [11]): 
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Proposition 4.1. A unitary n x n matrix U is permutation-symmetric if and only if U = al^^'^ + 
6J(") for a,b & C satisfying \a\ = 1 and \a + nb\ = 1. If U is moreover Hermitian, then U = 
±(jW _ |J(")). 

We see that only two permutation-symmetric Hermitian unitary matrices exist, both corre- 
sponding to d = f — 1- However, once the permutation symmetry is weakened to the modular 
permutation symmetry, there is much more freedom for d, as we shall see. 

In this section we will examine general properties of Hermitian unitary MPS matrices, in par- 
ticular necessary conditions of their existence, whereas sufficient conditions and concrete examples 
of such matrices will be presented in Section 5. 

Proposition 4.2. Let S e M„{d). Ifn>2, then d < f - 1. 

Proof. The diagonal entries of S arc +r and — r for r = --^^===. Since n > 2, at least two of 
them are equal, we may suppose without loss of generality that <Sii = 822- Moreover, we assume 
<Sii = +r; alternatively we would work with the equivalent matrix — <S. The unitarity of <S requires 
[SS*]i2 = 0, where 

n 

[<S<S*]^2 = <5ll«S21 +«Si2«S22 + ^^«Slj«S2j . 

i=3 

Let us denote Sjk = tc'"^'" for t = -^=5==.. Since S is Hermitian, it holds »S2i = »Si2, and thus 
the condition [SS* ]i2 = leads to 

n 

3=3 



hence we obtain 



therefore 



t 2 ^ 

J=3 



d=- <hn-2) = ^-l. 
t ~ 2^ ' 2 



□ 

Now we derive a relation between d and the signs of the diagonal entries of <S. 

Proposition 4.3. Let S G A4n{d), let p denote the number of its non-negative diagonal entries, 
and let m be the multiplicity of its eigenvalue 1 . Then 

2m-n={2p-n) . (6) 

Vm +n — 1 

Proof. Since <S e A4n{d), its diagonal entries are ±-^^7===^- According to the assumptions, 

^riS) = p-^^£— + {n-p) (- ^^4n-i ) = (2p-'^) vd4n-i - ^^'^ ^^^'^'^ ^ unitary 

and at the same time Hermitian, its eigenvalues are from the set {1, —1}, sec Observation 3.1. The 
multiplicity of 1 is m, the multiplicity of — 1 is n — m, hence Tr(iS) = m-l + (n — m) • (— 1) = 2m — n. 
Comparing these two expressions for Tr(<S) we obtain Equation (6). □ 

Notation 4.4. From now on to the end of the paper, the symbols m and p are reserved for the 
multiplicity of the eigenvalue 1 and the number of non-negative diagonal elements, respectively, of 
matrices <S S Mn{d). 



7 



Example 4.5. We demonstrate the use of formula (6) on the extremal values of d, namely d = 

(conference matrices) and d = f — 1- 

• Let d = Q. Then Equation (6) gives m—2n = 0. Consequently, complex Hermitian conference 
matrices exist only for even n. 

• Let d= f — 1. Then Equation (6), 2m — n = {2p — n) (l — has just three integer solution 
pairs: {m,p) = (1,0), {m,p) = (n— l,n), and {m,p) = §), the third one only for even n. 
These three solutions together with the parametrization of unitary matrices (2) can be used 
to an easy construction of all elements of A^„(^ — 1), cf. also [12]. 

Theorem 4.6. Let S G J^n{d) and •m,p have the usual meaning (see Notation 4-4)- Then 

• either p = m = ^ , 

• orp<m<^ orp>m> f; in this cased= |m- f | ^ (p_^)(~+^J ^- 

Proof. Any <S G M.n{d) satisfies Equation (6), which yields the following alternative: 

• 2m — n = 2p — n = 0, i.e., p = m = ^. 

• |2m— n| > |2p— n| > and at the same time = '^2p-n ■ ^® written equivalently 



as d = I'm — ^\ Y (p-rn)(p+m-n) ' where moreover p and m must satisfy p > m > f or 
p < m < f . 

□ 

Remark 4.7. Let S € Mn{d) for d ^ {|"^ " 1 1 \/ (p-m"(p+m ^ ' f < w < p < n|. Then, with 
regard to Theorem 4.6, m= p = which means in particular that n must be even. 
We finish the section with a remark on a matrices <S with m = 

Remark 4.8. Let S G M.n{d) for m = ^. Due to Theorem 3.2 (i), there exists a square matrix 
T G C"'™ such that 

^ +2 (/(™) +TT*)"^ 2 (/(™) +Tr*)"V \ 

I 2T* (J^™) + TT*)~"^ _/(™) _(_ 2T* (/(™) + TT*)^"^ T J 

Among all matrices of this type, those with normal T are particularly useful, because is such a 



case + 2T* (/('"^ + TT*) ^ T 



_/(™) +2 (/(™) -i-TT* 



F G 
G* -F 



, and consequently 

(7) 



where F = -/(") + 2 (/(™) + TT*) ^ and G = 2 (/(™) + TT*) V. Matrices S of type (7) are 
easier to be constructed, we will take advantage of this fact in the following section. 

5 Construction of Hermitian unitary MPS matrices 

Let us propose several ways how matrices M.n{d) can be constructed for certain values of d. 

First of all, for n = 2 there exists an <S e A42{d) for any d > 0, moreover <S can be always 



chosen real: <S = g^rpj ^ ) 
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From now on let n > 2. Now d is bounded from above by 2. — 1 (Prop. 4.2). With regard to 
this fact, we will structure our presentation according to the value of d, starting from the upper 
bound. The proposed matrix constructions mostly satisfy m = p = ^, and will be moreover based 
on the scheme! (7) from Remark 4.8; recall that any setting different from m = p = ^ would lead 
to a significant restriction on the admissible values of d, see Theorem 4.6. 



Case = f - 1 

There exists an S G A^n(§ — 1) for all n> 2, and S can be chosen real. 

According to [12], any S G M.nij ^ 1) is equivalent either to 7^"^ — ^J^"'^ or to the matrix 
obtained from (8) below by setting a = 0. 



Case de [f - 3, f - 1) 

There exists an S G M.n{d) for all even n gN. For example, set m = ^ and 

1 / (d + !)/('") - J(™) (e^" - + J('") \ 

^rf2 + „ _ 1 \^ (e-i« - 1)/M + JM -{d+ 1)7M + J^"*) J ' 

where a is chosen so that cos a = d + 2 — ^. 



Case d e (f - |, f - 3) 

// there exists an Hadamard matrix of order j + 1 or a symmetric conference matrix of order f + 1, 
then there exists an S G Mn{d) for all d G (f — |) f — 3) . 
To prove this statement, the following definition is needed. 

Definition 5.1. (i) Let H be an Hadamard matrix of order N having the form 



/ 1 



H 



1 



V 1 



1 



1 \ 



K 



H 



(9) 



The (A'' — 1) X {N — 1) matrix Kh is called the core of the Hadamard matrix H. 
(ii) Let C be a conference matrix of order N having the form 



/ 



c = 



\ 1 



1 



1 \ 



Kc 



(10) 



The {N — 1) X (N — 1) matrix Kc is called the core of the conference matrix C. 

It is easy to sec that if an Hadamard matrix of order N exists, then an Hadamard matrix of 
the form (9) exists. Similarly, if there is a conference matrix of order N , then there is a conference 
matrix having the form (10). 

Proposition 5.2. Let Kh he a core of a real Hadamard matrix of order f + 1. We define a matrix 
Gh of the size | x | as 



[Gh. 



jk 



ja[KHhk Q, chosen such that d = -2 - 



-(1 + cos^ a) . 
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Then 

o_ 1 / (rf +!)/("') -JM Gh \ 

V _(d+l)7M + jM; ^^^^ 

wif/i m = ^ satisfies S G Mn{d) for all d G [f ^ | - ^ l] • 

Proof. It suffices to prove tfiat 55* = J^"^. Witfi regard to (11), tliis condition is equivaient to 

1 f{d+lfI+[m-2{d+l)]J+GHG*H GhJ-JGh ^ = K") 

d^+n-l\ JG*h-G*hJ id+l)^I+[m-2id+l)]J + G*HGHj 

Since Kh is a core of an Hadamard matrix, it liolds: 

(a) every row and column of Kh contains tlie same number of I's and — I's, 

(b) the multiset {[KH]jk - [KH]ek \k = 1 . . . ,m} {j ^ £) equals {2, ■ ■ • . 2. -2. • • • , -2, 0, • • • , 0}. 

(m+l)/4 (m+l)/4 (m-l)/2 

From (a) it follows Gh-J-JGh = 0. From (b) we obtain [GhG*^^ = I2±i(e2i" + e-2i") + ™^ for 
all j + I. Since moreover {GhG*^]^^ = m, we have G^GJ^ = m/+ \^{(?'°' + e'^'") + ^\ ( J- 
/). These facts together with the relation between d and a lead to SS* = Finally, if a. runs 
over [0, tt], d defined by d = -2 + ^^(1 + cos^ a) runs over [f - |, f - l] • □ 

Remark 5.3. The construction (11) from Proposition 5.2 can be generalized. Let a symmetric 
(§, fc, A)-design exist and let A be its incidence matrix. We define K = ^(A + /) and a matrix 
G by Gj£ = e^"''^^^ for a satisfying d = -l + f-(fc-A)(l - cos 2a). Then the n x n matrix 
S constructed by analogy to (11) belongs to A^„((i). In this way one can obtain S E Ain{d) for 
d> —I + ^ — 2{k — A). Note, however, that for any symmetric {v, k, A)-design, the value fc — A is 
bounded from above by (see e.g. [13], Thm. 3.1.2), therefore always d > dmin = f ~ |) where 
dmin is the minimal value that has been obtained in Proposition 5.2. 

Proposition 5.4. Let Kc he a core of a symmetric conference matrix of order f + 1 and Gc be 
defined by [Gcjjk = e'"[^'^l^'=, where a is chosen such that d = + cos a + cos^ a. Then 



1 f {d+ 1)/('") - J^™) Gc 
+ V G*c -{d+ 1) J W + J W 



(12) 



n 3 1 n 1 

4 2 n-2' 2 ^ 



with TO = f satisfies S G Mn{d) for all d G 
Proof. Similar to the proof of Proposition 5.2. □ 

Case = f - 1 

There exists an S G Mn{d = f — |) for any even n. 

An S G Ain{d = f ~ |) can be constructed as follows. Let Kh be a core of a complex 
Hadamard matrix of order ^ + 1, i.e., Kh is of the size f x 2:. Note that there exists a complex 
Hadamard matrix of any order N, e.g. the one given by Hji~ = Q^'^'^(j-'^)0^-'^)/^ . Then 

with m = f satisfies S G Mn{j - §)• 
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Casede(t-|-^,t-i) 

// there exists a symmetric conference matrix of order f + 1, then there exists an S G A4n{d) for 
all de (^j - ^ - f - • See construction (12). 

Case d G [0, 1] 

If there is an Hermitian conference matrix C of order § ( equivalently: ^ (0) ^ 0^, then there 
exists an S £ A4n{d) for all d € [0, 1]. It is easy to check that this matrix can be constructed as 



+ C C 



e" 



Vrf2+n-l V - e--/(™) -(rf7('") + C) 



(14) 



where to = and a is chosen such that d = cos a. 



6 Real case (Symmetric orthogonal matrices) 



In this section we will focus on the matrices S € M.n{d) with the additional property that all their 
entries are real, i.e., on symmetric orthogonal matrices of the type 



n ■ 



±d 


±1 


±1 ••• 


±1 


±1 


±d 


±1 ••• 


±1 


±1 


±1 


±d ■■■ 


±1 


±1 


±1 


... ±1 


±d 



In what follows we will denote the real subset of A4n{d) by M^id), i.e. 
Ml{d) = 1^ e U(n) n M" 



5 is MPS A 



d A S =S^ 



r 



Elements of A^*(0) and A^^(l) represent (up to the factor -^=5===) symmetric Hadamard and 

symmetric conference matrices, respectively, of order n. For this reason, matrices S G Mn{d) with 
d G [0, f — 1] can be regarded as a straightforward generalization of the concept of symmetric 
Hadamard/conference matrices. A special subset of them, namely matrices <S € Al^(rf) with 
constant signs of the diagonal elements, have been studied in [2]. In this section we are interested 
in the case with general, mixed diagonal signs. 

Let us begin with examination of matrices of small orders. If n is small, it is an easy excercise to 
find admissible values of d using the orthogonality of the matrix rows. The results are summarized 
in the following Observation. 

Observation 6.1. It holds: 

• Ai^(d) is non-empty for all d E [0,oo); 

• Aif{d) is non-empty if and only if d ■ 

• (rf) is non-empty if and only if d=l. 

For dealing with <S e A4^{d) for a general n, let us introduce a notion of the standard form: 



1 . 

2' 
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Definition 6.2. We say that a matrix <S e M.^{d) is in the standard form if 





1 +d -1 •• 


• -1 




\ 


> 




-1 ■•• 








> 


1 


-1 


+d 














~d +1 •• 


• +1 










+1 ■•• 








I 




+1 


-d ) 





(15) 



n — p 



and p > J. For any S in the standard form, wc define Q ^ V d'^ + n — IS, and denote the blocks 
of g by QW,g("\g("^),Q(^^), where Q(^) is the left upper one of size p x p, i.e., 



1 



Vrf^ + n - 1 



If <S,<S e A1^((i), S ^ S and the matrix 5 is in the standard form, we say that <S is a standard 
form of S. 

Remark 6.3. Apparently, for any S G Ai^{n) there exists its standard form 5 ~ 5; on the other 
hand, such an <S is generally not unique. 

Lemma 6.4. Let S G Mn{d) be in the standard form (15) and let p > 3. 

(i) If there exist j, fc G {2, 3, • • • ,p}, j ^ k, such that Qj^j^ = +1, then n + 2d — 2 = (mod 4) 
and n — 6d — 6 > 0. 

(ii) If there exist j, k € {2, 3, • • • ,p}, j ^ k, such that Qj^^ = —1, then n — 2d — 2 = (mod 4). 
Proof. The first p rows of Q form the matrix 



( Q^'^ I Q^''^ ) 



/ +d 
-1 

V -1 



+d 



±1 ±1 



±1 \ 



since the rows of (g(^)|Q(")) are multiples of the rows of S, they are mutually orthogonal. 

For all j e {p+1, . . . ,n}, let us multiply the j-th column of {Q^^^\Q^^^^) by —Qij, which turns 
all the entries on the first row of Q^^^^ into —1. 

(i) Let = +1 for certain fc e {2, 3, • • • j ^ fc. In this case we apply the following 

two transpositions simultaneously to rows and columns of {Q'^^^\Q^^^^): 2 o j, 3 O /c. Note that 
this operation does not affect the orthogonality of the rows. As a result, the first three rows are 



d -1 -1 -1 1 -1 1 -1 1 -1 1 

-1 d (+1) +!■■■ + 1 +1--- + 1 -1 1 -1 1 

-1 (+1) d +1--- + 1 -1 1 +1--- + 1 -1 1 



(16) 
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They are orthogonal vectors from M^'", hence these four equations must be fulfilled: 

3 + £i + £2 + is + £4 = n (the vectors have n components) (17) 

-2d- 1 -4 -^2 +^3 +4 = (row 1 and row 2 are orthogonal) (18) 

—2d — 1 — ^1 + ^2 ^ ^,'5 + ^4 = (row 1 and row 3 are orthogonal) (19) 

1 + 2d + £1 - £2 - ^3 + ^4 = (row 2 and row 3 are orthogonal) (20) 

Wc Slim lip all the four equations to obtain 2 — 2d + ^£4 = n, and from (17)+(20) — (18)— (19) we 
get 6 + 6d + 4£i = n. Since £1 G No and £4 G No, it holds 

n + 2d - 2 = (mod 4) and n - 6d - 6 > . 

(ii) Let Qj-p = —1 for certain j, A; G {2, 3, • • • ,p}, j ^ k. Similarly as in the part (i), we apply 
the transpositions 2 ^ j , 3 ^ k simultaneously to rows and columns of (Q^^'lQ^"^) to rearrange 
the first three rows into the form 

d -1 -1 -1 1 -1 1 -1 1 -1 1 

-1 d (-1) +1--- + 1 +1--- + 1 -1 1 -1 1 

-1 (-1) d +1--- + 1 -1 1 +1--- + 1 -1 1 ^ ' 



In the same way as above, we obtain equations 6 — 6d + ^£4 = n and 2 + 2d + 4^i = n, hence 

n - 2d - 2 = (mod 4) . 

□ 

Theorem 6.5. If S e M^id) /or d < f - 1, then 

• n is even and n > 6, 

• d G No, 

• ^ + d is odd. 

Proof. Let S G A^*(d) for d < ^ — 1. Then, with regard to Observation 6.1, we have n > 5. We 
may assume without loss of generality (cf. Rem. 6.3) that <S is in the standard form. Therefore 
P > f (see Def. 6.2), hence p > 3, which allows us to use Lemma 6.4. Let Q, Q^^\Q^^^^ have the 
meaning introduced in Definition 6.2. We divide the explanation into three alternatives: 

• (The "positive" case.) Let us assume at first that Qj^^ = +1 for all j, fc G {2, 3, • • • ,p}, j ^ k. 
The orthogonality of the first two rows of <S gives the condition 

n 

-2d - (p - 2) + J2 QijQ-2j = . 
j=p+l 

However, since J2^=p+i QijQ'ij < n — p and at the same time it is assumed p > f , the 
condition cannot be satisfied for any d < f — 1- Consequently, the "positive" case is not 
possible. 
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• (The "mixed" case.) Let there exist j,k,j',k' G {2,3, ••• ,p}, j k, f 7^ k' such that 
'Sjfe' = +1 and QjVfc' = ~1- Then both statements (i) and (ii) of Lemma 6.4 apply, whence 
we get 

n + 2rf - 2 = (mod 4) and n - 2rf - 2 = (mod 4) . 

The first condition, being equivalent to ^ + = 1 (mod 2), means that ^ + is odd. Moreover, 
together with the second condition, it implies 2n — 4 = (mod 4) and 4(i = (mod 4), hence 
n is even and d is integer. 

• (The "negative" case.) Let finally Q^-^^ = —1 for all j,k S {2, 3,--- ,p}, j ^ k. Here we 
distinguish two situations: 

— If p = n, we have <S = l)/^") — J^"^). In this case, the orthogonality of 
the first two rows requires —2d + n— 2 = 0, hence d = f — 1, which contradicts our 

assumption d < f — 1- 

— If p < n, then the orthogonality of the 1st row and the (p + l)-st row of Q leads to the 
condition 

p n 

d • Qp+1,1 — Qp+i,j + Qi,p+i • i~d) + Qi,j=0- 

j=2 j=P+l 

Since Qp+1,1 = Qi,p+i, the terms with d cancel. The remaining condition is of the 
type 1 + 1 + -- - + 1— 1— I--- — 1 = for a certain £ S No, and thus can be satisfied 

e n-2-e 

only when n is even. Using this fact together with the relation n — 2d — 2 = (mod 4) 
from Lemma 6.4 (ii), we obtain d € No and ^ — d = 1 (mod 2), which is equivalent to 
f + d= 1 (mod 2). 

Finally, the inequality n > 6 follows from n > 5 and from the even parity of n, which has been 
proved above. □ 

Remark 6.6. It follows from Theorem 6.5: 

• If n is odd, necessarily d = ^ — 1. 

• If n EE 2 (mod 4), then d e {O, 2, 4, l}. 

• lfn = (mod 4), then {l, 3, 5, . . . , f - l}. 

Consequently, a real symmetric conference matrix (corresponding to d = 0) can exist only for 
n = 2 (mod 4). 

It turns out that for certain values of d, a more detailed description of <S can be found. Let us 
start with the following observation. 

Observation 6.7. The matrix J'^*-' has a simple eigenvalue k corresponding to the eigenvector 
w := (1, 1, ... , 1)-^, and the eigenvalue of multiplicity k — 1 corresponding to the eigenspace w-^. 

Proof. It holds J^'^^'w = kw, and since rank(j('^)) = 1, J^'^) has the eigenvalue with multiplicity 
dimker(jW) = fc-rank(jW) = fc- 1. □ 

Observation 6.7 will help us to characterize S G M^{d) for d exceeding f — 1: 
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Proposition 6.8. (i) Let S e M^{d) for d & (f - 1, f - l). Then p = ^ and there is a 
G e {-1,1}^'^' such that 



{d + - J(f) G 

V -id+ 1)/^^^ + J^"^ 



(22) 



The matrix G has these properties: G is normal, G commutes with J^^^ , and GG^ = {n — 2d — 
2)7(p) + (2d + 2-n/2)j(P). 

(a) On the other hand, if d and G fulfil the conditions above, then any S satisfying (22) belongs 

to Mn{d). 

Proof. We assume without loss of generality that S is in the standard form. Since n > 6 according 
to Theorem 6.5, it holds p > f > 3. Let Q^'\Q^"\Q<^'^^\Q<^^^'> have the meaning introduced in 
Definition 6.2. The proof will be carried out in five steps. 

Step 1. Since p > 3 and d > — 1, it immediately follows from Lemma 6.4 (i) that Q^^"* = 
(d+l)7(p) - J(P). 

Step 2. We prove that m > p. 

Since Q(i) = {d+ - J(P), it holds 

m = rank(5 + 7) > rank ^-^=^=Q(^) + /(f)^ = rank ((d + 1 + Vd^+n- l) I^p'> - J^^)) . 

Our goal is to show that the matrix L := (d + 1 + \/d'^ +n — l) 1^^^ — J^^^ is regular (equivalently, 
it docs not have the eigenvalue 0). 

Since the eigenvalues of 7*^^^ are and p according to Observation 6.7, the eigenvalues of L 
are d + 1 + \/d^ + n — 1 and d + 1 + Vd^ + n — 1 — p. The former one is trivially nonzero, thus it 
suffices to show d + 1 + Vd^ + n — 1 — p ^ 0. Wc proceed by contradiction: 

Let d + 1 + Vd^ + n — 1 — p = 0. Then the value £ := Vd^ + n — 1 — d = p— I — 2d satisfies 
f e N (because pGN and d e N, cf. Thm. 6.5), and also 

n— 1 n— In — 1 



£ = V d^ + n - 1 - d = —;==^=^ = < < 3 , 

Vd^ + n - 1 + d p-l-f-1 

because n > 6. Consequently, £= \ ox £ = 2. 

• Case £= \ implies \/d^ + n — 1 — d = 1, hence d = f — 1, which contradicts the assumption 

|_l<d<|_l. 

• Case £ = 2 implies Vd^ + n — 1 — d = 2, hence d = However, since n is even (see 
Thm. 6.5), such d is not an integer, and thus is not admissible. 

Therefore L is regular, hence rank(L) = p, thus indeed m = rank(L) > p. 

Step 3. Since p > f and at the same time m > p due to Step 2, Theorem 4.6 implies 
m = p = |. 

S'tep ^. We show that = -(d + 1)J("-p) + J("-p). 

Since p = § and S is in the standard form, it is obvious that the matrix S' := —PSP~^ 



i^p^ 

for P = I ^(p) Q 1 , which takes the form S' = 



satisfies 



<S' ~ 5 and is in the standard form as well. Therefore, with regard to Step 1, it holds —Q^^^^ = 
(d + l)/(p) - J(P). 

Step 5. The result of Step 1 together with the hermiticity of S imply (22). From the unitarity 
of S, SS* = we immediately obtain the properties GG^ = G^G, GJ^p^ = J^p^G and 

GG'^ = (n - 2d - 2)/(p) + (2d + 2 - n/2)j(p\ And vice versa, if G fulfils these conditions and 
d^ (t ~ 1' f ~ 1)' the matrix (22) satisfies SS* =1. □ 
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Proposition 6.8 allows us to formulate the necessary and sufficient condition of the existence of 

S e Ml{d) for all d > f - 1: 

Theorem 6.9. (i) There is no S € Ml{d) /or d e (f - 1, f - |) . 
(a) An S e A4n{d) for d G [f~i)f~l) exists if and only if the value 

is an integer and there exists a symmetric (f , A;, X) -design for k = j — ^ and A = ^{d—q+1). 

Proof, (i) According to Proposition 6.8, an <S e M^{d) for d > ^ — 1 exists if and only if there 
exists a normal G € {— 1, {p= satisfying 

Gj{p) = j(p)g , GG'^ = {n-2d- 2)1^^^ + (2rf + 2 - n/2) J(p) . 

With regard to Observation 6.7, the matrix GG^ has a simple eigenvalue ^ + — l){2d + 2 — ^) 
and a corresponding eigenvector w = {1,1, . . . ,iy . Since GG^ is a nonnegative matrix, necessarily 
§ + (f - l)(2d + 2 - I) > 0, hence we obtain the condition 

d>^--^--^. 
- 4 2 n-2 

Finally, wc know from Theorem 6.5 that d, is integer, n is even and n > 6, for this reason the last 
condition can be equivalently written as d > — |. Henee (i) is proved. 

In the proof of (ii), we start from the following claim: 
Claim. If M is normal, then v is an eigenvector of M with a simple eigenvalue a if and only if v 
is an eigenvector of MM* with a simple eigenvalue |crp. (The proof is straightforward using the 
eigendecomposition of M.) 

We have found in part (i) that w = {1,1, . . . , 1)^ is an eigenvector of GG^ corresponding to a simple 
eigenvalue ^ + {^ — l){2d+2 — ^) > 0. Therefore, due to Claim, it is at the same time an eigenvector 

of G corresponding to a simple eigenvalue of modulus q := y^f + (f — l) {"^d + 2 — ^) . Since 
both w and G are real and their entries are integers ±1, the eigenvalue n must be real and integer, 
hence q G Nq and /i = ±q. The actual sign of /i is irrelevant, because we can always turn G in (22) 
into — G by multiplying the rows ^ + 1, . . . ,n ol S and the columns + 1, . . . , n of 5 by — 1. Let 
us assume for definiteness /x = -q- We define A=^{G + J(p)); then A e {0, 1}^'^ and 

^ « ' ^ ' . ' 

'= fe-A ^ 

It follows (see (1)) that A is the incidence matrix of a symmetric (f , A;, A) -design. □ 

Remark 6.10. After finishing this work we found out about paper [3] in which orthogonal (not 
necessarily symmetric) MPS matrices with constant integral diagonal have been studied. It follows 
from there that our Theorem 6.5 is valid even when the symmetry of the matrix is weakened to 
the non-skew-symmetry. 

Let us finish the section by a series of remarks on construction of matrices <S e A4^{d). 
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Notes on constructions of <S e -M^id) 

• AnS e A^*(f - 1) exists for all n>2. For example S = - f J^"\ see also Section 5. 

• An S e — 3) exists for any even n. To obtain <S, set a = tt in (8). 

• AnS € Ml{d) ford€ [f - §, f - 3) exists iff q := y^f + (f - l) (2d + 2- f) e No and 

there exists a symmetric (f , fc, A) -design for k — j — ^ and A = ^{d— q+1). (Cf. Thm. 6.9.) 

The construction follows from Theorem 6.9. If A is the incidence matrix of a symmetric 
(i'?~2'i~2 + 5) -design, we set G = 2^ - j(t) and construct <S according to (22). 

• In particular, for = ^ — | we obtain: 

Observation 6.11. There exists an S G — |) i/ and only if there exists a real 

Hadamard matrix of order f + 1 • 

Proof We apply Theorem 6.9. If = f - |, then g = 1 e Nq and A: = f - i = ^ (f + l) - 1, 

A=|-— | = -j(f + l) — 1. Therefore the existence of an iS e Af*(^ — |) is equivalent to the 
existence of an {N — 1, — 1, jN — l)-design where N = f + 1, which is further equivalent 
to the existence of an Hadamard matrix of order N (see [14], Lemma 1.9.3). □ 

But we could prove Proposition 6.11 in a direct way as well: 

The validity of implication <= is confirmed by the construction (13), therefore it suffices 

to prove =>. Let S € A^*(^ — |). Due to Proposition 6.8, S satisfies (22) for a certain 
G e {—1,1}^'^. In the same way as in the proof of Proposition 6.8, we can show that 
w = {1,1, . . . , 1)^ is an eigenvector of G corresponding to a simple eigenvalue fj, = ±1. We 
define 





1 ••• 1 \ 


1 






G 


1 


J 



and using the properties of G derived in Proposition 6.8, we show that HH* = (^ + 1)7, 
thus H is an Hadamard matrix. 

Remark. It follows from Theorem 6.5 that an<S € — |) exists only if n = 6 (mod 8). 

Observation 6.11 implies that if the Hadamard conjecture is true, then an S G A^*(5 ~ 5) 

exists if and only if n = 6 (mod 8). 

• //d e (f - 1, f - I), no 5 e Ml{d) can be constructed. (Thm. 6.9.) 

• // there is a symmetric conference matrix C of order ^, then there exists an S G A4^{1). To 

obtain S, set a = in (14). 

• The existence of iS S A^*(0) is trivially equivalent to the existence of a symmetric conference 
matrix of order n. 

• For constructions of <S € M^{d) with p = n, we refer to [2], where real symmetric MPS 
matrices with constant diagonal were studied and certain methods of their construction have 
been proposed. 
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7 An application: Quantum graphs 



Let us briefly explain in what context matrices from the sets A4n{d) emerge in quantum mechanics 

on graphs. 

Consider a metric graph, i.e., a set of vertices and a set of edges, the edges connect the vertices, 
each edge has a given length. Let us suppose that the graph is of microscopic size and that there 
is a particle, for example an electron, having certain energy and moving along the graph edges. 
As the size of the system is very small, the behaviour of the particle is governed by the laws of 
quantum mechanics. In particular, its position cannot be exactly determined, one can only flnd 
the probability density of its occurrence in a given point x of the graph, which is given as |\E'(a;)p, 
where \1/ is the wave function of the particle. The function depends on the topology of the 
graph, on the lengths of the edges, on the particle energy E and on physical characteristics of the 
vertices (junctions). The physical characteristics of each junction are expressed by the scattering 
matrix S that has the following properties: 

• 5 is a complex n x n matrix, where n is the vertex degree. 

• Let the edges coupled at the junction be numbered by 1, . . . , n. If the quantum particle comes 
in the junction from the j-th line, then it is scattered into all lines 1, . . . ,n (including the 
j-th line itself) with the probabilities . . . , |iS„jp. In other words, the squared moduli 
of the entries of <S correspond to the scattering probabilities at the junction. 

• S is always unitary (this property may be viewed as the quantum version of Kirchhoff 's law, 

see [15]). 

• S generally depends on energy (where S{E) can be uniquely calculated from <S(1)). 

• <S is energy-independent if and only if <S is Hermitian. 

Now it is obvious what role the Hermitian unitary MPS matrices play in this theory. Consider 
a junction of degree n. If its physical characteristics are described by a Hermitian unitary MPS 
matrix <S G Mn{d), then the particle is transmitted from any edge to any other edge with equal 
probabilities, also the reflection probabilities are the same at all edges, and furthermore, the 
probabilities are independent of the particle energy. The parameter d squared represents the ratio 
of the reflection probability to the probability that the particle is transmitted to any chosen edge 
different from the incoming one. 

It is also noteworthy that real scattering matrices S, examined in Section 6, correspond to 
junctions with the additional physical property of time reversibility. 

Let us add that the graph with a particle is a model of realistic physical systems where a 
particle moves along thin nano-sized wires made for example of semiconductors. That is why 
the study of the existence of Hermitian unitary MPS matrices is important - the existence of 
an 5 € Ain{d) determines whether or not it is possible to physically construct (manufacture) an 
equally-transmitting junction with the given scattering ratio d^. 
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